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Generalized Parton Distributions - theoretical review - 
D. Miiller a 

a Institut fiir Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, 
Germany 

In this talk an introduction to generalized parton distributions is given. Recent devel- 
opments are shortly reviewed, including non-perturbative calculations, phenomenological 
aspects and evaluation of higher order perturbative and power corrections. 

1. Introduction 

Generalized parton distributions (GPDs) [ HJ El El an d their crossed version, i.e., the 
generalized distribution amplitudes (GDAs) [HJE] appear in the perturbative description 
of certain hard exclusive processes, e.g., in the leptoproduction of photon and mesons 
or in the production of hadron pairs in two photon annihilation. Although the first 
systematic study of GPDs and GDAs was done for more than one decade [ [Q, their 
physical significance has been widely realized in connection with the proton spin puzzle. 
Namely, these distributions incorporate, besides other non-perturbative information that 
is not encoded in forward parton densities or form factors, also gravitational form factors 
from which the quark orbital angular momentum fraction contributing to the nucleon spin 
can be read off [|2|. This perception induced then intensive theoretical and experimental 
studies, in which it has been fully realized that GPDs and GDAs are a new concept to 
study the structure of hadrons and nuclei, contain a link between exclusive and inclusive 
processes and open a new window for the exploration of non-perturbative QCD. In the 
following a mini review about these developments is given, for comprehensive ones see 
Refs. [000!. 

In Sect. 121 the basic properties of GPDs and their partonic interpretation are presented. 
In Sect. E]GPD ansatze are shortly discussed and is then devoted to the non-perturbative 
evaluation. The exclusive processes, which allow to gain access to these distributions, are 
listed in Sect. HI In Sect. results from the evaluation of radiative and power suppressed 
corrections beyond leading order (LO) are reported and finally, the conclusions are given. 

2. Features of Generalized parton distributions 

GPDs are defined as Fourier transform of light-ray operators, sandwiched between the 
initial and final hadronic states. There is a whole compendium of GPDs: corresponding 
to the species of hadrons, the initial and final states can have different quantum numbers 
(transition GPDs), and even we might replace the hadrons by nuclei (nucleus GPDs). 
Specified the initial and final states, GPDs are classified with respect to the twist of the 
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operators and the spin content of fields. At leading twist-two level three different types 
of quark and gluon GPDs can be defined (here the gauge link is omitted): 
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with P + = n ■ (Pi + P 2 ), V- = n* ■ V, n 2 = (n*) 2 — 0, n ■ n* — 1. In the first (vector) and 
second (axial- vector) entry the in- and outgoing partons have the same helicities, where 
its sum or difference of left- and right-handed partons is taken, respectively. For the 
third entry, called transversity, a helicity flip appears. GPDs depend on the momentum 
fraction x, conjugated to the light-cone distance 2k, the longitudinal momentum fraction 
£ = (Pi — P 2 ) + / (Pi + P 2 ) + in the t-channel, the momentum transfer A 2 = t = (P 2 — Pi) 2 , 
and the renormalization scale fi 2 . The latter is induced by the renormalization prescription 
of the operators, which is part of the GPD definition. To deal with the polarization of the 
hadronic states, one might introduces a form factor decomposition [El GUI- For instance, 
for the nucleon GPD Dirac and Pauli-like form factors appear in the vector case 

\ V = Z7(P 2 , S 2 ) 1+ U(P U SjH^x, £, A 2 ) + U(P 2 , S 2 )^±^P(Pi, S 1 )E i (x, £, A 2 ) , (3) 

where i = u,d, s, ■ ■ ■ ,G. Definitions (HJ-© imply the basic properties: 

• GPDs reduce in the forward limit A — > to parton densities [ [TJ El 121 Hj , e.g., 

qi(x,fi 2 ) = Jm5,(x,(,A 2 ,|i 2 ), and hm E^x, f, A 2 , /i 2 ) ^ 0, (4) 

while the helicity flip contribution Ei decouples. 

• The /i 2 -dependence is governed by linear evolution equations [[I], which can be 
derived from the renormalization group equation of the light-ray operators [ ITTl ll2j. 

• Hermiticity together with time reversal invariance leads to a definite symmetry with 
respect to the skewness parameter £, e.g., Pj(x,£) = Hi(x, — £). 

• The Mellin moments of GPDs are expectation values of local twist-two operators: 
J dx x n \ v {x, e, A 2 , Q 2 ) = n« • • -n^(P 2 , 5 2 |S^ 7w) iS w • • -W^ V q \Pi, Si) . (5) 

Lorentz covariance induces that they are polynomials in £, which are even in (JSJ. 

Furthermore, GPDs are constrained in the region x > |£| by the positivity of the norm 
in the Hilbert space of states. The most general form of such positivity bounds [ I13| I14j. 
known so far, are given as an infinite set of constraints [I15j. Since GPDs are implicitly 
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Figure 1. Partonic interpretation and support property (right) of GPDs. 



scheme dependent, positivity bounds are a request on both the GPD model and scheme. 
Note also that GPDs can be represented as overlap of light-cone wave functions [ ■ 

Let us give a partonic interpretation of GPDs. A generic quark GPD, e.g., in the vector 
case, q(x, £, A 2 ) = Q(x, £, A 2 ) — Q(— x, £, A 2 ) is decomposed in its quark Q and anti-quark 
Q part. In the central (exclusive or ER-BL) region — £ < x < £, Q(x,£,A 2 ) might be 
interpreted as probability amplitude to have a meson like configuration inside the hadron, 
while the outer (inclusive or DGLAP) region £ < x < 1 can be viewed as probability 
amplitude for emission and absorbing a quark with momentum fraction X\P\ = j^Pi 
and x<iP\ = f^f-fi; respectively. Lorentz invariance ties this both regions, which can be 
read off from the representation 1 that ensures polynomiality (valid for £ > 0): 

Q(x, e, A 2 ) = 9 (-£ < x < 1) la, (|, i A 2 ) - 9 (£ < x < 1) i ro -i A 2 ) . (6) 



Obviously, in the central region the GPD is given by w from which the outer region, 

determined by the antisymmetric function w f|, j \ — w {— |, — |J, can be restored, see 
Fig. The uniqueness of this continuation procedure was shown in connection with the 
support extension of evolution kernels [[!]. 

GPDs simultaneously possess a longitudinal and transversal momenta dependence and 
so they encode the three dimensional distribution of partons in the considered hadron or 
nucleus [ I17j. Indeed, a partonic density interpretation holds in the infinite momentum 
frame as long as £ = [I18j. see also Refs. [I19| I20j. In this kinematics the central region, 
i.e., the parton number violating contributions, drops out. Thus, it could be shown that 
in the impact parameter space within the infinite momentum frame 

K(x,b ± ) = J*^e*^H i {x,£ = 0,-Al) (7) 

has the probabilistic interpretation to find a parton species % in dependence on the mo- 
mentum fraction x and the relative distance b± from the proton center. It is noted that 
an interpretation of the three dimensional Fourier transform of GPDs in the rest frame 
has been suggested within the concept of phase space (Wigner) distributions [ |2H 122] ■ 

^q. where zu(x,y, A 2 ) = J +v dw d(x,x — my, A 2 ) might be derived by means of a partonic Fock 
state decomposition and the so-called a-representation for Feynman diagrams. 
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3. Parametrization, model and lattice calculations of GPDs 

GPDs are mostly unknown functions with a complex variable dependence, which must 
satisfy the basic properties. Because of the lack of further knowledge, GPD ansatze, 
needed for the estimation of cross sections and asymmetries, are constructed so far by 
simplicity and intuition [Ej. Such an ansatz, given at the input scale Q 2 ,, is based on 
the assumption of (x, £) and A 2 factorization and is now widely used in phenomenology: 

Hi{x, £, A 2 , Ql) = F,(A 2 )^(x, , hi(x, = f 1 dy 5(x - y - f z)D i (y, z) . (8) 

J-l J-l+\y\ 

Here Fj(A 2 ) are partonic form factors, partially fixed by sum rules. The reduced GPD 
hi(x,£) is given in terms of a factorized ansatz for the so-called double distribution 
Di(y,z) = qi(yQl)U(\z\/(l - \y\))/(l - \y\) with the parton density qi(y,Ql) and an 
unknown profile function II [ El IE] ■ This ansatz might be refined by a regge-like ansatz 
for the parton densities at small x: qi(x) ~ x~ ai — > x ~ ai ~ a i A [|H]. To repair an artifact 
in the relation among GPD and DD [[22] a so-called D-term is added in (J5|). intuitively 
understood as meson exchange contribution in the t-channel [l2*lj. 

In the region £ < x for small x it also has been proposed to equate the reduced 
GPDs with the parton densities at the input scale and so skewness is purely generated 
by evolution We note that in this region, where gluon and sea quark contributions 

are dominating, the parameterization of the hard-scattering amplitude convoluted with a 
GPD drastically simplifies. Here, the skewness effect at the input scale effectively enters 
the normalization [I2E]- A further proposal for the GPD parametrization is based on their 
representation as an infinite sum of t-channel exchange contributions [I27j. 

The GPD parameterizations, proposed so far, do not simultaneously satisfy the needs: 

• Basic properties must be automatically fulfilled. 

• Flexible parametrization of the degrees of freedom that are left. 

• Simple numerical treatment of evolution and convolutions. 

These requests can be mostly satisfied within a generalization of the Mellin representa- 
tion for parton densities by adopting the concept of complex angular momentum to the 
conformal spin expansion of GPDs [12*5]. 

The poor knowledge about GPDs might be improved by non-perturbative model calcu- 
lations, done at a low fi 2 scale. The resulting GPDs can then be evolved to the scale that 
is relevant for their phenomenological use. However, the scheme dependence of GPDs 
leads to uncertainties in this matching procedure, which should be considered as part of 
the model. Because of limited space, all the efforts cannot be reviewed here in detail, 
however, at least it should be mentioned a few of them, namely, calculations within the 
MIT bag model [ 29, chiral quark soliton model [ED], constituent quark model [EI], and 
Bethe-Salpeter and light-cone wave function approaches [E2]- 

More recently, the first few Mellin moments (jSJ) for proton GPDs have been measured 
on the lattice [EH]- Especially, the quark orbital angular momentum fraction of the proton 
spin, the second moment of a certain GPD combination, could be extracted. Also the 
correlation of the A 2 dependence with the order n of moments leads to a valuable insight 




Figure 2. Factorization of the meson leptoproduction amplitude (left) and kinematics. 



in the nucleon GPDs. The slope of the A 2 dependence decreases with increasing n and 
so the transversal size of the proton shrinks with increasing x, which has been argued in 
Ref . [ • From this one concludes that the factorized ansatz (JHJ) is oversimplified. It 
might be also possible to rich from lattice results a better understanding of the skewness 
dependence. So far the common GPD ansatze leads always to an enhancement of the 
(reduced) GPDs, compared to parton densities, at the crossing point x = £. For the 
phenomenology it is highly desired to have a more rigorous understanding of this issue. 

4. Hard exclusive reactions to leading order 

GPDs appear in the perturbative description of the hard photon or meson leptoproduc- 
tion: l^N — > l^BX, where B is the final baryon state and X stands for the photon or the 
observed meson. The virtuality — Q 2 of the intermediate photon must be sufficiently large. 
This implies that the hadronic scattering amplitude factorizes into a hard-scattering one 
convoluted with GPDs and eventually a meson distribution amplitude (DA), see Fig. El 
[ EJ EH HHH EH] ■ Usually, one defines this factorization in such a way that the collinear 
divergencies are removed within a (modified) minimal subtraction scheme. Obviously, 
this procedure is ambiguous and induces, e.g., the factorization scale dependence. If the 
same factorization scheme is applied for all processes, GPDs and DAs are universal in the 
sense that they are process independent, however, they depend on the scheme and of the 
order at which the perturbation theory is truncated. Thus, GPDs serve at the first place 
as a tool that connect physical observables measured in different processes. Concerning 
their probabilistic interpretation, one must conclude that the number of measured par- 
tons depends on the scheme conventions (even in untruncated perturbation theory). This 
problem, appearing also for parton densities, might be resolved by choosing a reference 
process, e.g., the leptoproduction of photons, in which all radiative corrections (order by 
order) are absorbed in the GPD definitions. 

Let us first consider the leptoproduction of a meson. So far the quantum numbers 
allow such a process and the virtual photon is longitudinally polarized, the process is 
for sufficient large Q 2 perturbatively described as shown in Fig. El The initial and final 
hadron state might have different quantum numbers and so it is applicable for several 
processes, namely, the production of neutral and charged vector and pseudo scalar [ 
E2H mesons and or even for exotic states like the pentaquark or hybrid mesons [I41j. 
Besides the cross section also the transversal target spin asymmetry is described within 
the perturbative framework, which might be justified at a lower scale Q 2 as for the cross 
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Figure 3. Deeply virtual Compton scattering and Bethe-Heitler bremsstrahlungs process 
(left) and factorization of the Compton amplitude to leading order (right). 

section itself [ 021 S3] • To LO the convolution of GPDs and meson DA are separated 



where qff and 4> A j refers to the (transition) GPD and meson DA, respectively. Qf,Qf are 
the partonic charge factors, corresponding to the flavor content of the observed meson. 
Hence, the produced meson serves as a flavor filter and so a flavor decomposition for GPDs 
can be reached, e.g., by measuring the processes e~p — > e~pir and e~p — > e~pr]. The 
amplitude Al scales with 1/Q, which leads to a scaling law da/dA 2 oc l/<2 6 for the cross 
section. However, this canonical scaling is logarithmically modified by evolution effects, 
mainly arising from GPDs. 

Several of these processes have been measured or a planned in fixed target (HERMES, 
JLAB, COMPASS) and collider (HI, ZEUS) experiments, for a review see the contribu- 
tions of F. Sabatie and L. Favart in this proceedings. I only like to give a short attention 
to the electroproduction of neutral (longitudinally polarized) vector mesons, measured in 
the small x^ ~ 2£ region up to rather large Q 2 by the HI and ZEUS collaborations. Here 
the amplitude (jHJ) drastically simplifies, since now it is dominated by gluon exchange, and 
thus such experiments are an ideal testing ground of the perturbative framework and the 
study of the gluonic GPD. Let me remind that under the assumption of SU(3) symmetric 
meson DAs the ratio of cross sections should be: a p o : : = 9 : 1 : 2, which is in fair 
agreement with experimental data (plotted with respect to the scale Q 2 + M v ). 

Another class of hard exclusive processes, in which GPDs are accessible, are those in 
which the struck parton is purely probed within the electromagnetic interaction via the 
absorbtion and emission of photons [H]: photo leptoproduction [113122], see Fig. El lepton 
pair photoproduction [HSJ and lepton pair leptoproduction [ESI- In the first two processes 
the photon virtuality is space- and time-like respectively, while in the third one both 
photons are virtual. The latter one is experimentally most challenging, however, it is the 
only one in which a deconvolution of the GPDs is possible. In all these processes there are 
two interfering subprocesses, the hard virtual Compton scattering and the Bethe-Heitler 
(BH) bremsstrahlungs process, given in terms of the electromagnetic form factors. For 
the case that at least one photon has a large virtuality, the Compton scattering process 
factorizes, shown in Fig. |H]to LO. 

Hence, the cross section has a rich azimuthal angular dependence and the interference 
term is linear in the GPDs. This offers the possibility to unreveal the GPDs by measur- 
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ing asymmetries that are dominated by the interference term: charge, single beam and 
target spin asymmetries, while double spin asymmetries require the subtraction of the 
BH contribution. In collider experiments it is even possible to extract the deeply virtual 
Compton scattering (DVCS) cross section, since the BH amplitude is sufficiently small. 
A detailed compendium of asymmetries and their relations to GPDs is given in Ref. [ 26 . 
The first measurements of single beam spin and charge asymmetries and cross sections 
are compatible with the oversimplified GPD ansatz (jHJ). 

It was also argued that the amplitudes of the photon and meson photoproduction at 
large s and —t factorize to LO in terms of the inverse GPD moments j\dxq(x,£,t)/x 
[ El SHI ESj- For these processes also a perturbative treatment is used in which all 
valence partons are resolved via hard gluon exchanges [ E7| • In this conjecture here 
only the struck parton is resolved while all spectators are contained in the GPD. This 
is somehow analogous to the mechanism proposed by Feynman for the description of 
elastic form factors. Unfortunately, no factorization theorem could be established, since 
power corrections are uncontrolled. This problem has been attempted to resolve by the 
assumption that the virtuality of the partons in the initial and final state is always small 

Note that the perturbative framework is also applicable for the photon [ED] and meson 
[ I5T] leptoproduction off nuclei. This opens a new window for the study of nucleus binding 
effects, especially, for the deuteron [ [TU1 E2] • 

5. Beyond leading order predictions 

To gain insight in the validity of the perturbative framework and to improve it, it is 
necessary to calculate higher order perturbative and power suppressed corrections or at 
least to estimate them. A large amount of work has been done in this direction during 
the last few years. 

The factorization theorems, derived to leading power accuracy, state that the perturba- 
tive corrections are systematically calculable in this approximation. The next-to-leading 
order (NLO) corrections for the hard-scattering amplitude of both Compton scattering 
processes [ EH1 EEJ an d f° r leptoproduction of mesons have been completed [ HH1 EI] . 
While for the former process the perturbative corrections, which depends on the GPD 
ansatz, turn out to be moderate, the size of NLO corrections for the latter is rather large. 
However, these corrections partly cancel in the transversal proton spin asymmetry [ I43j . 

Employing conformal consistency predictions [EH]' the evolution kernels have been 
evaluated to NLO [EH] and are implemented in numerical codes [Ell- Note that conformal 
symmetry in the minimal subtraction scheme is broken in a subtle manner, which can be 
removed by a finite renormalization, providing the conformal subtraction (CS) scheme [ 
EH]- In such a scheme the next-to-next-to-leading order (NNLO) corrections to both the 
DVCS hard-scattering amplitude and evolution kernels can be borrowed from the known 
results in deep inelastic scattering, for a first discussion see Ref. [155]. 

Whether power suppressed contributions can be calculated within perturbative QCD 
depends on the process in question. For hard exclusive leptoproduction the exchange 
of a transversal polarized photon yields 1/Q suppressed contributions. However, they 
are affected by non-integrable singularities, which appear to LO in a s and are induced by 
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large size quark-antiquark configurations of the produced meson [Hj. In the case of photon 
leptoproduction 1/Q suppressed contributions are calculable to LO [ ISO] and NLO [EI] 
accuracy. In fact, they complete the azimuthal angular dependence of the cross section [ 
Interesting to remark, that within an appropriate definition of observables they do 
not interfere with the twist-two prediction, appearing at leading power. 

Twist-two predictions will be affected by 1/Q 2 power suppressed contributions appear- 
ing at twist-four level. Such contributions have been perturbatively calculated to LO by 
neglecting multi-particle operators [E21- Unfortunately, this approximation suffers from 
an ambiguity in the choice of the multi-particle operator basis, which shows up in the 
violation of current conservation. This problem is not resolved so far. Upper bounds 
of the power corrections to the hard electroproduction have been estimated within the 
renormalon approach and found to be large [I21j. 

The fact that perturbative and non-perturbative corrections for leptoproduction of 
mesons are much larger as for photons indicates that the onset of the scaling region 
for the former process will be at a higher scale as for the latter one. For photon leptopro- 
duction one would expect that, as in the case of deep inelastic scattering, a scale of few 
GeV 2 is sufficient to apply perturbative QCD. 

6. Conclusions 

In the last decade the concept of generalized parton distributions and distribution am- 
plitudes has been enormously developed: a partonic interpretation has been given in 
depth, factorization theorems has been derived, predictions for numerous processes has 
been worked out, radiative corrections are calculated to NLO accuracy, including the evo- 
lution equations, first measurements on the lattice has been performed, and experimental 
accessibility has been demonstrated in pioneering experiments. 

GPDs and also GDAs can be explored within measurements of hard exclusive, however, 
inelastic processes in fixed target and collider experiments, where a deconvolution with 
respect to the momentum fraction is practically impossible. Only the leptoproduction 
of a lepton pairs allows to scan the central region of GPDs. These functions are hy- 
brids that incorporate parton densities, form factors, and distribution amplitudes. Thus, 
they form a link between exclusive and inclusive processes and encode information about 
non-perturbative QCD, which cannot be obtained from inclusive or elastic exclusive pro- 
cesses. GPDs/GDAs are a new tool that allows to study the structure of hadrons and 
nuclei from a new perspective. So for instance, their knowledge would provide the quark 
orbital angular momentum contributing to the hadron spin and the three dimensional 
distribution of partons inside hadrons and nuclei. Moreover, they can serve for the study 
of flavor and chiral symmetry breaking or nucleus binding effects, to name a few. As it 
has been stressed, GPDs are scheme dependent quantities and, thus, their probabilistic 
interpretation should be done within an appropriate scheme convention. 

The solution of several open problems will require a large effort. The most challenging 
issue for the theory is whether factorization is applicable for the photon virtuality reached 
in present experiments. Here higher order calculations beyond the NLO accuracy might 
give some insight. Also the growing amount of experimental data will help to answer 
this question. A further issue concerns the appropriate and realistic parametrization of 
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GPDs. From the first photon leptoproduction measurements it becomes also clear that 
high precision data are needed to distinguish between GPD models. Certainly, in the last 
decade a huge, however, first step has been taken in reaching a deeper understanding of 
the hadron and nuclei structure within the concept of GPDs/GDAs. 
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